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We show that the electrical impedance of a small-capacitance Josephson junction includes besides 
the capacitive term —i/uiCB also an inductive term iluLb- Similar to the known Bloch capacitance 
Cs(g), the Bloch inductance Ls(g) also depends periodically on the quasicharge q, and its maximum 
value achieved at q = e (mod2e) always exceeds the value of the Josephson inductance of this junction 
Lj(if) at fixed tp — 0. The effect of the Bloch inductance on the dynamics of a single junction and 
a one-dimensional array is described. 
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In the recent years the Josephson tunnel junctions 
with small capacitance have been extensively studied. 
The interest in these junctions and in circuits with these 
junctions can be explained by the remarkable quantum 
behavior of these macroscopic systems. For example, 
the small-capacitance junctions can transfer individual 
Cooper pairs and generate Bloch oscillations 0, 0, 0|- 
These junctions can also exhibit a long-time quantum co- 
herence 0, 0] • Important potential applications of these 
effects include the fundamental standard of current 6] 
and quantum computing circuits [7|. 

The small-capacitance Josephson junctions are charac- 
terized by the finite ratio A = Ej/E c of the Josephson 
coupling energy Ej = (<E> /27r)/ c (here I c is the critical 
current and <f>o = h/2e is the flux quantum) and the 
charging energy E c — e 2 /2C (C is the junction capaci- 
tance). The remarkable feature of these junctions is their 
nonlinear differential capacitance Cb related to the local 
curvature of the zero Bloch energy band Q, , 
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where E (q) is the ground state energy of the junction. 
Eq periodically (the period is equal to 2e) depends on 
the quasicharge q — J I(t')dt' , which is a good variable 
driven by the current source I (see Fig. la). Variable q 
is analog to the quasimomentum of electron moving in 
the periodic potential of a crystal lattice. Similar to the 
Josephson inductance Lj(ip), which depends on the clas- 
sical phase ip in the large junctions (A — > oo), capacitance 
Cb can also take on both positive and negative values. 
The property of variation of the Bloch capacitance in 
a large range was turned to advantage by Averin and 
Bruder || who suggested a switchable capacitive cou- 
pling between Josephson charge qubits. Incorporating 
the nonlinear Bloch capacitance into a tank circuit made 
it possible to realize a sensitive electrometer 111 and effi- 
cient readout of the Cooper-pair-box qubit [lTJ . 

Behavior of these circuits is usually described by the 
junction model comprising solely the Bloch capacitance 
Eq. QJ. Actually, this approach is applicable for the 
static case or for a sufficiently slow variation of the 



quasicharge q. In this paper we show that the small- 
capacitance junction model should generally include an 
inductance term (see Fig. lb). Such a term is propor- 
tional to the second time-derivative q — I, so the corre- 
sponding inductance, which we will call the "Bloch in- 
ductance," plays the role of the mass for quasicharge q. 
We evaluate the effect of the Bloch inductance and show 
that it can be essential even for the single-band motion 
of the system. 

The circuit consisting of a small-capacitance Josephson 
junction biased by the ideal source of the classical current 
I(t) is described by the Hamiltonian H = H + Hr 
where the junction's Hamiltonian and the term due to 
the source are respectively equal to 



and ffj = -|V(t). ( 2 ) 



The operator of the charge Q is conjugate to the phase 
operator <p and is equal to Q = ^£ -Q- in (^-representation. 
The quasiparticle tunneling is negFected for the sake of 
clarity. The eigenenergies E n (k) and the Bloch eigen- 
states \k,n) of the base Hamiltonial Ho are periodic func- 
tions of the continuous quasi- wavenumber k — q/2e or 
the quasicharge q, with the period equal to 1 and 2e, 
respectively; the band index n = 0, 1, 2, ... [ll| . 

The operator of the phase is given in this basis by the 
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FIG. 1: (a) Electric diagram of the current-biased small- 
capacitance Josephson junction and (b) the equivalent circuit 
of this junction for a small ac signal. Connected in series 
the Bloch capacitance Cb and Bloch inductance Lb both de- 
pend on the instant value of the quasicharge q(t) = J 1 1(t')dt' 
which is set by the current source. 



2 



matrix 



elements 00 02 



V>1& = ig^5(k - k')5 nn , + V l n '8{k - k')(l - <W)- (3) 

Since Hamiltonian H is diagonal with respect to the vari- 
able k, the density matrix elements p^, take the form 
(Tnn>5(k — k'). After taking trace over this variable, the 
equation of motion p = (l/ih)[H, p] is reduced to 
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Here w fc m = [i? n (fc) — E n i(k)\/h is the interband angular 
frequency, // is the instant Bloch frequency fj — I(t)/2e, 
and the quasi-wavenumber time dependence is governed 
by h, 



k(t) = k(0) + / fl(j)dT. 
JO 



(5) 



The operator of voltage, V = (&fJ2n)<p = (i/2e)[H,(p], 
is given by the matrix elements |12l Il3| 



1 dE n 
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The diagonal term Voo = V(k) gives the voltage value 
in the ground state (see the plots for a wide range of 
parameter A in Fig. 2b of Ref. 0|)> while its derivative 
gives the reverse capacitance value Cg 1 Eq. Q shown in 
Fig. 2a. The observable value of the voltage is equal to 



(V)=Tr{Vp}= V nn>a n >n, 

n,n'—0 



(7) 



so the non-zero off-diagonal terms of the voltage operator 
can also, in principle, contribute to (V). 

We solve Eq. J2J taking into account the two lowest 
bands, i.e. n,n' — and 1. We also assume that the ex- 
citation to the state n = 1 is small, so (Too i5 F while 
fiii l cr oiU <J iol ^ 1- This assumption is adequate if 
both the instant Bloch frequency fi(t) and the maxi- 
mum rate of its alteration fi(t) are sufficiently small, i.e. 
the dimensionless parameters a = max(//)/f2fc <C 1 and 
ol\ = max(//)/f2^ <§; 1, where the transition frequency 



is denoted as fife = a 
equations of motion 



,10 



In this case, the solution of the 



r = o"oo - &i\ = 2if I ((p° k 1 a w - (p{°a i), (8) 



&oi = v*o = -i^fecroi - ifiy°k r ( 9 ) 



reads 



<Toi(t)=o-Io(*)=e- </ » n *" , * / «roi(0) 
_ ie -iJ*n h ,dt' f e i/ t 'n fc „ ( jt"^oi /sJ(t /)^/ ) (n) 



where we introduce the quasi-wavenumbers k' = k(t') 
and k" = kit"). Keeping the first-order contribution 
with respect to a and a%, we finally arrive at 
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The first term describes the time-evolution (oscillations) 
of the density matrix element due to Hamiltonian Hq. 
The amplitude A of these high-frequency oscillations de- 
pends on the initial conditions, i.e. 



A = M0) + fM//(0) + i 



n 



(13) 



where t^g 1 = v f(t = 0) and f2 = &k(t = 0). For the 
system initially prepared in the ground state, a nn i (0) = 
<5n,o^o,n', and at "smooth" switching-on of the current 
I(t), i.e. 1(0) = and 1(0) = 0, the amplitude A = 0. 
The last term on the right-hand side of Eq. (|12J) is, how- 
ever, essential, yielding the dependence of eroi on fj = k. 
The observable value of the voltage Eq. (JTJ is equal to 



(V) = Voodoo + 2K(y 10 a i 
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and can therefore be presented in the form 



(V) — V(k) + L B (k) 



dl 
dt ' 



(15) 



Here we introduce the Bloch inductance 



L B (k) 



E 



J -\vf\ 2 Lj^ (16) 



where Ljq is the Josephson inductance of a similar classi- 
cal junction (i.e. in the absence of charging effects, E c — > 
0), taken at zero argument, Lj((p = 0) = <I>o/(27r/ c ). One 
can see that in contrast to the phase-dependent Joseph- 
son inductance Lj(ip), the Bloch inductance L B (k) is a 
positive periodic function of the quasi-wavenumber (qua- 
sicharge). Physically, L B still characterizes the kinetic 
properties of supercurrent which is now the operator [l3| 
whose value depends on the system state, i.e. on the good 
variable k, or, equivalently, a very wide wave packet over 
phase ip in (^-representation. 

Let us consider the behavior of the Bloch inductance at 
small and large A. In these limiting cases, the analytical 
expressions for the matrix elements ip® 1 can be taken from 
Ref. 0. For A <C 1, inductance L B increases resonantly 
in the vicinity of the degeneracy points, i.e. 



coo = 1, o"n = 0, 



(10) 
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FIG. 2: Normalized reverse Bloch capacitance (a) and Bloch 
inductance (b) in the zero band (ground state) versus qua- 
sicharge for several values of the energy ratio A = Ej/E c . 



at k 

T max 
B 



(1 + A£)/2 (mod 1). The peak value is equal to 
32A~ 2 L,/o 3> Lj . Outside the resonance region, 
Lb drops dramatically approaching at k « (mod 1) the 
lowest level of Uf a = 0.5A 3 L J0 < L J0 . 

At large Josephson couplings, A ;§> 1, both the 
interband frequency ilk « uj p = (8EjE c y/ 2 /Ti — 
(S/Xy^Ej/h (ujp is the Josephson plasma frequency) 
and the matrix elements l^ 1 « |<Poscl = (2/A) 1 / 4 (equal 
to those of the equivalent harmonic oscillator) do not 
depend on k. The Bloch inductance takes on almost 
the constant value, Lb ~ Ljq. Taking into account 
the fact that corrections to ip®l c are exponentially small 
(oc exp [— (8A) 1 / 2 ] ) and accounting only for corrections to 
the oscillator's energy eigenvalues due to anharmonicity 
of the Josephson potential, the expression for Lb reads 



L B (k)= [l + (2A)- 1 / 2 



L 



jo ■ 



(18) 



For several intermediate values of parameter A, the induc- 
tance Lb is presented by the plots in Fig. 2b, obtained 
by numerical methods. 

If the dc part of the quasi-wavenumber k = q/2e is 
fixed and the small ac part is driven by current I(t) = 
Iie luJt , the junction voltage according to Eq. I|15|) is equal 
to V = V(q) + Z(u>)Iie ZUJt . The electrical impedance of 
the junction Z(uj) = i[loLb — (wCs)" 1 ] is purely reac- 



tive (see Fig. lb). So, at a sufficiently high irradiation 
frequency uj, the inductive term can be appreciable. 

At small Josephson couplings, A <C 1, the maximum 
effect of Lb can be expected to be near the degeneracy 
points Eq. lfT7|) . However, even in these points the induc- 
tive resistance luL b approaches the capacitive resistance 
— (wCb) -1 (the positive value due to the negative value 
of Cb = —AC/4) only at ui = ilk = Ej/h, i.e. at the 
frequency of the interband transition. At such a high 
frequency our approach is not applicable. At lower fre- 
quencies, the effect of Lb is rather small and can be pre- 
sented as a frequency-dependent correction to the reverse 
Bloch capacitance, 



r- 1 - 

Cb XC 



~xc + M 



(19) 



Such renormalization of the Bloch capacitance can, 
in principle, be detected in experiment on the 
high-frequency B loch-capacitance-based readout of the 
Cooper pair box qubit [lfj. 

For A » 1, as follows from the harmonic shape of the 
eigenenergy E (q) (see Eqs. (A3)-(A4) of Ref. 0), the re- 
verse Bloch capacitance is small, 



Cg\k) = C^l cos(27rfc), C B l = bB(X) C 



B{X) = A 3/4 exp [-(8A) 1 / 2 



«T 



(20) 



(21) 



where the numeric factor b = 7r 3 / 2 2 11 / 4 « 37. The max- 
imum (resonance) frequency above which the inductor 
dominates, is 



= [L B (0)C B0 }- 1/2 = (bB)^ 2 uj p « n k . 



(22) 



Therefore, at least in the range of the frequencies ujq < 
uj -C ilk) the effect of inductance Lb should be signifi- 
cant. 

The question arises whether the role of the Bloch in- 
ductance is essential in the dynamics of an autonomous 
Josephson junction. For evaluating this effect we intro- 
duce a small dissipation by adding to our model Eq. 
a high-ohmic resistor, R ^> Rq = h/4e 2 « 6.45 kfi, 
connected either parallel to the current source in Fig. 1 
0,0 or, equivalently, in series with a voltage source Vq, 
as shown in the inset in Fig. 3. In this case the sys- 
tem is described by a wave packet over k (see Eqs. (71) 
and (75) of Ref. 0). At low temperature, k B T < 
mm{hflk, [Eo(e) — Eo(0)]}, this packet with the center 
k, is sufficiently narrow, Ak -C 1. Corresponding inter- 
band damping terms in the equation of motion Eq. @ 
lead to an exponential decay of the fast-oscillating off- 
diagonal term cx A in Eq. i|12[l . Finally, the equation of 
motion for q takes the form 



L B (q)q + Rq + V(q)=V , 



(23) 
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or, in dimensionless units, 



(24) 



Here we introduced 8 = 2irk, r — w c t, lu c — V c /eR, 
v = Vo/Vc, V c = max[V(q)] and the unity-amplitude pe- 
riodic function g(6) = V(e8/ir)/V c . Parameter (3l = 
LjoV c /eR 2 and the normalized Bloch inductance is 
£(8) = L B (e9/-K)/Lj . For large A, the values of these 
parameters are V c ~ n^ 1 eC] 3 Q, uj c » (RCbo) -1 , 9 (9) ~ 
sin 9, i{8) w 1 and 
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Rq_ 
R 



« 1, (25) 



and Eq. Q24|l takes the form of the equation, describ- 
ing a driven, damped pendulum or a resistively shunted 
Josephson junction 15]. Finite values of corresponding 
McCumber-Stewart parameter, /3j, > 1, characterizing 
the effect of the inductive term on the dc I-V curve 
(IVC), can therefore not be achieved. 

However, the situation changes dramatically in the 
case of an array of N junctions connected in series in 
the same network (see the central part of the inset in 
Fig. 3). The individual junctions of this array are decou- 
pled and described by Eqs. (|23|1 - (|25|) with replacement 
L B (q) -» NL B (q), V(q) -» NV(q), V c -» NV C and 
(3l — > fi\ — N 2 Pl. Note that the maximum resonance 
frequency ujq Eq. H24[l does not change. At sufficiently 
large N 3> 1, the effective McCumber-Stewart parame- 
ter P* L can be made of the order of 1, even for unavoidably 
large values of R. As a result, the IVC of such an array 
can exhibit characteristic hysteresis (see Fig. 3). 

In experiment, due to the small stray capacitance Cq of 
each island of the array to ground (see the rightmost part 
in the inset in Fig. 3), the number of decoupled junctions 
iVo is limited by the size of the Cooper-pair soliton, i.e. 
N = (C B0 /C ) 1/2 [Hill- At A > 1, this size increases 
dramatically, so that it can result in hysteretic IVC. Sim- 
ilar effect has been observed repeatedly by KTH group in 
long arrays of small Al junctions (see, e.g., H3)j although 
these arrays were voltage-biased, so the quasiparticle tun- 
nelling played an important role in those samples. 

The concept of the Bloch inductance allows the 
sine-Gordon equation, describing static Cooper-pair 2e- 
solitons in the small-capacitance junction array |l6| . to 
be generalized to the non-stationary case. Each junc- 
tion included in the elementary cell of the array is pre- 
sented as shown in Fig. 1. Due to the emerging inertial 
term L B q, the corresponding sine-Gordon equation for 
the quasicharge describes a persistent motion of the 2e- 
solitons. 

In conclusion, we introduced the Bloch inductance and 
described its effect on the quasicharge dynamics in cir- 
cuits with small-capacitance Josephson junctions. The 
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FIG. 3: I-V characteristics obtained in the model including 
the quasicharge-dependent Bloch inductance Eqs. 123H - 124|I 
for several values of parameter /3£ . The curves corresponding 
to the finite values of f3* L show appreciable hysteresis. The 
inset shows possible realizations of large values of /3£. 



effect should manifest itself by an effective renormaliza- 
tion of the Bloch capacitance in the vicinity of the degen- 
eracy points of the Cooper pair boxes and by a hysteretic 
behavior of the IVCs of long one-dimensional arrays. The 
Bloch inductance of the Josephson junction included in 
the variable electrostatic transformer [j| may at A ^> 1 
result in an inter-qubit coupling of resonator-mediated 
type. Finally, similar to the Josephson inductance the 
Bloch inductance may be applied for constructing the 
dispersive QND-type readout of Josephson qubits |l9| . 
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